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The decays J/tp —¥ pp and J/tp — ► nn have been investigated with a sample of 225.2 million J/tp 
events collected with the BESIII detector at the BEPCII e + e~ collider. The branching fractions 
are determined to be B(J/ip -» pp) = (2.112 ± 0.004 ± 0.031) x 10~ 3 and B(J/ip -» nn) = (2.07 ± 
0.01 ± 0.17) x 10~ 3 . Distributions of the angle 8 between the proton or anti-neutron and the beam 
direction are well described by the form 1 + a cos 2 8, and we find ct — 0.595 ± 0.012 ± 0.015 for 
J/tp — > pp and ct — 0.50 ± 0.04 ± 0.21 for J/tp —¥ nn. Our branching- fraction results suggest a large 
phase angle between the strong and electromagnetic amplitudes describing the J/tp — s> NN decay. 
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I. INTRODUCTION 



The J/ip meson is interpreted as a bound state of a charmed quark and a charmed antiquark (cc). The decay process 
J/ip — > NN (N = p or n) is an octet-baryon-pair decay mode, and should be a good laboratory for testing perturbative 
QCD (pQCD) because the three gluons in the OZI-violating strong decay correspond to the three qq pairs that form 
the final-state nucleons. The ratio of the branching fractions for the pp and nn final states provides information about 
the phase angle between the strong and the electromagnetic (EM) amplitudes governing the decay [IHU . Because the 
initial-state isospin is 0, the strong-decay amplitudes for the pp and nn final states must be equal. The J/ip — > PP 
and J/ip nn EM decays are expected to have amplitudes that are of about the same magnitude, but with opposite 
signs, like the magnetic moments (as discussed in Section fVII[) . Because the EM decays of J/ip to pp and nn behave 
the same as non-resonant production of those final states, the magnitude of the EM decay amplitude of J/ip can be 
estimated from the cross section for continuum production e + e~ — ¥ pp. If the strong and EM amplitudes are almost 
real, and therefore in phase, as predicted by pQCD [U-IH, then interference would lead to a branching fraction for 
J/ip — s- nn about one-half as large as that for J/ip — >■ pp. Conversely, if the strong and EM amplitudes are orthogonal, 
then the strong decay dominates and the branching fractions are expected to be equal. In previous experiments, 
JJip — > PP has been measured with good precision, while J/ip —> nn has been measured with quite a large uncertainty 
[a LZ|- They appear to be equal within errors, at odds with the pQCD expectation. 

The angular distribution for J/ip — > NN can be written as a function of the angle 6 between the nucleon or 
antinuclcon direction and the beam as follows: 



= A 1 +acos 2 , 

dcosO 

where A is an overall normalization. These angular distributions reflect details of the baryon structure and have the 
potential to distinguish among different theoretical models [l|-[j| . 

In this paper, we report new studies of the process J/ip — >• NN made with the BESIII detector at the BEPCII 
electron-positron storage ring With the world's largest sample of J/ip decays, we obtain improved measurements 

for the J/ip — >■ pp and J /ip — > nfi branching fractions and angular distributions. 



II. BEPCII AND BESIII 



BEPCII is a two-ring e + e~ collider designed for a peak luminosity of 10 33 cm _2 s _1 at a beam current of 0.93 A. 
The cylindrical core of the BESIII detector consists of a hclium-gas-based drift chamber (MDC) for charged-particlc 
tracking and particle identification by dE/dx, a plastic scintillator time-of- flight system (TOF) for additional particle 
identification, and a 6240-crystal CsI(Tl) Electromagnetic Calorimeter (EMC) for electron identification and photon 
detection. These components are all enclosed in a superconducting solcnoidal magnet providing a 1.0- T magnetic field. 
The solenoid is supported by an octagonal flux-return yoke with resistive-plate-counter muon detector modules (MU) 
interleaved with steel. The geometrical acceptance for charged tracks and photons is 93% of 47r, and the resolutions 
for charged-track momentum and photon energy at 1 GeV are 0.5% and 2.5%, respectively. More details on the 
features and capabilities of BESIII are provided in Ref. Q . 



III. DATA SAMPLE 



Our data sample consists of 225.2 million e + e~ — » J/ip events collected during 2009. The estimated uncertainty in 
the number of events is ±1.3% [l0|. A GEANT4-based [III [III detector simulation is used to produce Monte Carlo (MC) 
samples for signal and background processes that are generated with specialized models that have been packaged and 
customized for BESIII [13( • EvtGen |14| is used to study phase-space signal events for J/ip — > pp and for exclusive 
backgrounds in J/ip decays. BABAYAGA [l5| is used to generate Bhabha and 77 events as possible EM backgrounds. A 
large inclusive sam ple (200 million events) is used to simulate hadronic background processes. The J/ip resonance is 
generated by KKMC [17|. Known J/ip decay modes are generated with EvtGen, using branching fractions set to world- 
average values . The remaining J / ip decay modes are generated by LUNDCHARM [l3| , which is based on JETSET fl6j 
and tuned for the charm-energy region. The decays J/ip — > pp and J/ip — > nn are excluded from this sample. 
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IV. GENERAL EVENT SELECTION 



Charged tracks in BESIII are reconstructed from MDC hits. To optimize the momentum measurement, we select 
tracks in the polar angle range | cos#| < 0.93 and require that they pass within ±10 cm of the interaction point in 
the beam direction and within ±1 cm in the plane perpendicular to the beam. 

Electromagnetic showers are reconstructed by clustering EMC crystal energies. Efficiency and energy resolution 
are improved by including energy deposits in nearby TOF counters. Showers used in selecting photons and in 7r° 
reconstruction must satisfy fiducial and shower-quality requirements. Showers in the barrel region (| cos#| < 0.8) must 
have a minimum energy of 25 MeV, while those in the endcaps (0.86 < | cos$| < 0.92) must have at least 50 MeV. 
Showers in the region between the barrel and endcap are poorly reconstructed and are excluded. To eliminate showers 
from charged particles, a photon must be separated by at least 10° from any charged track. EMC timing requirements 
suppress electronic noise and energy deposits unrelated to the event. 



V. ANALYSIS OF J/ip ->■ pp 
A. Event Selection 



Events with exactly two good charged tracks in the polar angle range | cos#| < 0.8 are selected. We exclude the two 
endcap regions to reduce systematic uncertainties in tracking and particle identification. By using a loose particle- 
identification requirement for the positive track (probability of the p hypothesis greater than the probabilities for the 
7r + and K + hypotheses), and by requiring no particle identification for the negative track, the efficiency is maximized 
and the systematic uncertainty is minimized. A vertex fit is performed to the two selected tracks to improve the 
momentum resolution, and the angle between the p and p is required to be greater than 178°. Finally, for both tracks, 
the measured momentum magnitude must be within 30 MeV/c (~ 3<r) of the expected value of 1.232 GcV/c. Fig. [T] 
shows comparisons between data and MC for the angle between the p and p and for their momenta. 

This selection results in a signal of N = 314651±561 candidate events. Fig.Q](b) and (c) show the p and p momentum 
distributions for these events, along with the expected distributions for a pure MC J/tp — > pp signal. Backgrounds 
overall are very small, and appear to be negligible in the accepted p and p momentum range. Three independent 
procedures are used to estimate this background: inclusive J/ip MC, exclusive MC of potential background processes 
(Bhabha events and J/ip decays to e + e _ , /^ + ^~, K + K~ , 7pp, 7r°p|>, and 7?y c with rj c — > pp), and a sideband technique. 
The estimates range from 0.02% to 0.2% of the signal. We apply no subtraction and take the largest of the estimates 
(sideband) as a systematic uncertainty in the final result. The raw distribution of cos (9 for the protons in the selected 
signal events is given in Fig. [5J 




FIG. 1. Comparisons between data (points) and MC (histograms) for properties of the p and p tracks for selected J/if) — > pp 
signal events: (a) angle between the p and p, (b) p momentum, and (c) p momentum. 
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FIG. 2. Angular distribution of the selected J/ip — > pp candidates. 



B. Efficiency Correction 



To measure the J/ip —> PP branching fraction and angular distribution, it is necessary to correct for the selection 
efficiency, which is dominated by track reconstruction and selection and by particle-identification efficiency. We use 
signal MC to obtain the efficiency, but use data to correct for imperfections in the simulation, thereby reducing 
the systematic uncertainty in the correction. We measure differences between the data and MC separately for the 
efficiencies of tracking and particle identification, leaving the other selection cuts in place or tightening them for 
cleaner selection. The correlation between the corrections in the tracking and particle-identification efficiencies has 
been shown in MC studies to be small, so we combine them into a single correction function that is applied to the 
MC-determined efficiency. Because the tracking and TOF response depend on the track direction, the efficiency 
correction is determined in bins of cos 9. 

We divide the full angular range (| cos#| < 0.8) into 16 equal bins and for each bin compute the efficiency for the 
successful reconstruction of the p or p track as follows: 

N 2 

etrk ~ N 1+ N 2 > 

where N\ (N 2 ) is the number of J/tp —> pp events with 1 (2) good charged track(s) detected. For Ni, we require 
only one good charged track which is identified as a p or p. Note that in this case, unlike the J/ij; —¥ pp selection, we 
can apply particle identification to the p selection to improve purity, since any inconsistency between data and MC 
would cancel in the efficiency. Fig. [3] shows the data/MC comparison for the tracking efficiencies and the computed 
correction factor efrk a / e trk f° r eacn cos # bin for p and p. 
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FIG. 3. (a) The proton tracking efficiency for data (points) and MC (circles), and (b) the correction ttrk 1 l e u-k i ( c ) an d (d) 
show the same for antiprotons. 

We can similarly measure the particle-identification efficiency for the p in each cos# bin, considering only J/ip — > pp 
events in which there are two good charged tracks, with the negatively-charged track identified as an antiproton. We 
define the efficiency as follows: 



where N p is the number of selected events in which the proton has been successfully identified and Np is the number 
of events without the proton identified. To select a more pure sample we tighten the selection on the p and p momenta 
to be within 20 McV/c (~ 2a) of the expected value. Figure 2] shows the data/MC comparison for the proton particle 



identification efficiency and the resulting correction factor e^?j a ' /tpld ■ 

In each cos0 bin, the corrected MC-determined efficiency to be applied to data is computed with the following 
formula: 



e MC x 



-data -data -data 



ptrk 
c ptrk 



ppid 
ppid 



ptrk 

Taw 

ptrk 



To diminish the effect of bin-to-bin scatter due to statistical fluctuations, we fit the corrected efficiency as a function 
of cos 9 with a fifth-order polynomial, as shown in Fig. [S] 



C. Angular Distribution and Branching Fraction 



We fit the measured angular distribution of the proton from J/tp —> pp to the function A(l + a cos 2 8)e(cos8), 
where A is the overall normalization and e(cos#) is the corrected MC-determined efficiency function (Sect. fVBj) . The 
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FIG. 4. (a) p particle-identification efficiency for data (points) and MC (circles), and (b) the computed efficiency correction 
factor eggV#£. 
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FIG. 5. J/ip — > pp selection efficiency as a function of cos# (a) before correction, and (b) after correction. The line shows the 
smoothed efficiency obtained by fitting the data to a fifth-order polynomial. 




FIG. 6. The points represent the measured distribution of cos for the p in J/ tp —> pp candidate events, with error bars that are 
the quadratic sums of the statistical and efficiency uncertainties. The line represents the fit of the distribution to the functional 
form given in the text, and is used to determine the normalization and the angular-distribution parameter a. 
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angular distribution and the fit are shown in Fig. [6] The \ 2 f° r the fit is 16, with 14 degrees of freedom, and the 
value determined for the angular-distribution parameter is a = 0.595 ± 0.012, where the error is statistical only. 

The raw yield of J/ip —> pp events obtained by counting protons in the angular range cos 6* = [—0.8,0.8] is 
7V(— 0.8, 0.8) = 314651 ± 561. The efficiency-corrected yield obtained by fitting the cos 8 distribution over this range 
is N cor (— 0.8, 0.8) = 357786 ± 638. The fitted value of a is used to determine the total number of J/ip — > nn events 
in the full angular range of cos 9 = [—1.0, 1.0] as follows: 

f^Jl + acos 2 9)d9 

N cor (-1.0, 1.0) = N cor (-0.8, 0.8) x — = 475567 ± 848. 

J° l 8 (l + a cos 2 9)d9 

Combining this final yield with the number of J/ip events in our sample ((2.252±0.029) x 10 8 ), we find the branching 
fraction to be 

B{J/ip -> pp) = (2.112 ± 0.004) x 10~ 3 , 

where the error is statistical only. 

D. Systematic Errors and Results 

To determine the uncertainty in the efficiency correction, we use toy MC experiments to obtain distributions in the 
branching fraction and a that reflect the statistical errors of the bin-by-bin efficiency values. We perform this study 
by varying each bin randomly according to a normal distribution for each MC experiment, redoing the polynomial fit 
and then remeasuring the efficiency-corrected yield. The results have normal distributions and they are fitted with 
Gaussian functions to estimate the associated uncertainties in the branching fraction and a, which are found to be 
4.69 x 10~ 6 and 0.011, respectively. 

The full magnitude of the pp momentum sideband background estimate (0.2%) is taken to be the uncertainty in 
the branching fraction due to the background correction. We fit the sideband-subtracted angular distribution and 
determine a new value for the angular parameter a, taking the change relative to the standard result (0.004) as the 
systematic error. 

To estimate the systematic error due to the detector angular resolution, we perform a study with the signal MC. 
The "true" generated proton cos 9 is fitted before and after smearing with a MC-derived angular resolution function. 
The differences in the fitted a values (0.010) and in the branching fractions (4 x 10 -6 ) are taken as the systematic 
uncertainties from this source. 

The branching fraction also incurs two systematic uncertainties that do not affect a. A small systematic uncertainty 
enters due to the correction for the | cos#| < 0.8 requirement, which depends on the determined value of a and its 
error. The dominant uncertainty in the branching fraction is due to the estimated 1.3% error in the number of J/ip 
events in our sample [l(| ■ 

To study the effect from continuum production, we write the total cross section a p p as 



V FF s — m z + iml tot 

where E p and S are the EM and strong amplitudes of J/ip — > pp and </> is the phase angle between them. erSg n *' is 
the pp cross section contributed by the continuum under the J/ip peak. These values are taken from the calculation 
in Sect. IVII1 The difference with and without cr^™*- is assigned as the systematic error. 

Finally, we change a by ±1 a (includes systematic error) and reevaluate the branching fraction to estimate the 
systematic error in the branching fraction measurement. 

Table [J provides a summary of all identified sources of systematic uncertainty, which are assumed to be uncorrelated, 
and their quadrature sum. The final results for our J/ip — > pp measurements are as follows: 

a = 0.595 ±0.012 ±0.015, and 



B(J/ip pp) = (2.112 ± 0.004 ± 0.031) x 10~ 3 . 
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The branching fraction measurement is consistent with the previous world average [6j and improves the overall 
precision by about a factor of 2.5. The value of a is also consistent with previous experiments (Table [TTJ) and is 
improved significantly. 



TABLE I. Systematic errors for J/ip — > pp. 



Sources 


Effect on a 


Effect on B(10 -3 ) 


Efficiency Correction 


O.Off 


0.005 


Background 


0.004 


0.002 


cos 9 Resolution 


O.OfO 


0.004 


a Value 




0.004 


Number of J/ip 




0.026 


continuum 




0.015 


Total 


0.015 


0.031 



TABLE II. Previous measurements of a in J/ip — > pp. 



Collaboration 


a 


Markl [18] 
Mark2 [19] 
Mark3 [20J 
DASP [21] 
DM2 [22] 
BESII [23J 


1.45 ±0.56 
0.61 ±0.23 
0.58 ±0.14 
1.70 ± 1.70 
0.62 ±0.11 
0.676 ± 0.055 



VI. ANALYSIS OF J/ip -+ nn 
A. Event Selection 



We search for J/ip — > nn candidates by selecting events that have no good charged tracks originating in the 
interaction region. The antineutron annihilation "star" in the EMC provides a signature for these events that is much 
more identifiable than the hadronic shower produced by a neutron. We therefore first select events with showers 
characteristic of n interactions, and then search in these events for energy deposited by n hadronic interactions on 
the opposite side of the detector. 

The most energetic shower in the event is assigned to be the n candidate and is required to have an energy in the 
range 0.6 — 2.0 GeV. To optimize the discrimination against backgrounds, we apply a fiducial cut of | cos 0\ < 0.8 to the 
n candidate. This ensures that the n energy is fully contained in the EMC for most signal events. To suppress photon 
backgrounds, we impose a requirement on the second moment of the candidate shower, defined as S = TuE^r? /Ej-Ej, 
where Ei is the energy deposited in the i th crystal of the shower and Ti is the distance from the center of that crystal to 
the center of the shower. To be accepted, the n candidate must satisfy S > 20 cm 2 . To further exploit the distinctive 
n shower topology, we require the number of EMC hits in a 50° cone around the n candidate shower direction to be 
greater than 40. 

Events with accepted n candidates are searched for EMC showers on the opposite side of the detector that are 
consistent with being the neutron in a J/ip — > nn decay. The energy of this shower must be between 0.06 and 
0.6 GeV, a range found to be characteristic of the EMC neutron response in MC studies. If multiple showers arc 
present, the one that is most back-to-back with respect to the n candidate is selected. To further suppress backgrounds 
from all-neutral J/ip decays, continuum production and EM processes, we require E extra = 0, where E extra is the 
total deposited energy in the EMC, excluding that of the n shower and any additional energy in the 50° cone. 

The expected signal for J/ip — > nn is an enhancement near 180° in the distribution of the angle between the n 
shower and the direction of the n. The distributions of this angle and of the cosine of the polar angle of the n shower 
(cos 9) for selected candidates are shown in Fig.[7J The enhancement near 180° in the distribution of the angle between 
the n and h constitutes the J ftp — > nh signal. Since there is nonncgligiblc background, the number of J/ip — > nn 
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events must be determined by fitting. Distributions of the angle between the n and h are constructed in bins of cos 9 
and fitted with signal and background functions. 

Data-driven methods are used to determine the efficiency and signal shapes for J/ijj — > nh. We select J ftp — > p(n)Tt~ 
and charge- conjugate (c.c.) events in data to obtain n and h samples to evaluate the selection efficiency. We use the p 
and p in J /ip — > pp events that have been selected using information from just the MDC to get unbiased information on 
the shape and efficiency of the n and h response in the EMC, since antiproton and antineutron hadronic interactions 
are similar. 

Generic J/ip MC is used to assess the background. Fig. [7] (a) shows that there is no peaking in the distribution 
of the angle between the n and h for this background. We also consider possible exclusive background channels: 
J/ip — > ir°nh, J/tp — > jnn, e + e~ — > 77, J/ip — > J/ip — > J/ip — > pp, and J/tp — > 7?y c (n c — > nh). None 

of these potential background sources exhibits peaking in the distribution of the angle between the n and h. 
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FIG. 7. Distributions for J/tp — > nn candidate events summed over all cos 6 bins (points) and background from inclusive MC 
(solid lines): (a) angle between the n and n, and (b) cos# for the n shower. 



B. Efficiency Determination 

We use specially chosen event samples from data to determine the efficiencies for each requirement in the J /ip — > nn 
selection. The overall efficiency is then computed bin-by-bin in cos 6 as the product of these components and is applied 
as a correction in obtaining the angular distribution and branching fraction. 

We select J/ip — > prm~ events to study the efficiency of the fi selection. Events with exactly two good charged 
tracks identified as p and ir~ arc selected. Information from the TOF detector and dE/dx information from the MDC 
are combined to do the particle identification. The p and tt~ are required to have a missing mass within 30 McV 
of the nominal fi mass. The missing momentum of the p and tt~ is required to be in the range 1.1 — 1.2 GeV/c to 
ensure a sample that is as similar as possible to the n in J/ip — > nh (momentum 1.232 GeV/c). The number of events 
passing the above selection gives N exp , the expected h yield. The number of h candidates selected from these events 
(criteria defined in Sect. IVI A[) that match the missing momenta of the accompanying p and ir~ within 10° gives the 
observed yield 7V obs . The efficiency for n selection is ef ta = N° bs /N exp . 
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To validate this procedure and ensure consistency between the n in J/ip — > pfm~ and that in the signal process 
J/tp — > nn, we select higher-purity n candidates in J/ip — > nh (J/tp — > pn-K~) with a stringent cut of 177° on the 
angle between the n and n. (For J/ip — > pfnr~ the cut is on the angle between the n and the missing momentum of 
the p and 7r~.) Comparisons of the selection variables (energy deposit in EMC, number of EMC hits in a 50° cone 
about the n shower, and the shower second moment) for these two n samples are shown in Fig. [8] Each distribution 
is plotted after the cuts on the other variables have been imposed. There is good agreement, verifying that the n in 
J/ip — > p(h)n~ is a good match to the n in the signal process J/tp — > nn, and that this process provides a reliable 
efficiency correction. 

We apply the same technique to calculate the efficiency for selecting the neutron (e^ ata ), in this case using a 
sample of J/tp — > p(n)7r + selected from data. A comparison of the distribution of the EMC energy for neutrons from 
J/ip — > p(n)ir + with that from J/ip — > nn is shown in Fig. |9] In this case the momentum difference between the two 
n samples results in a greater difference in the EMC energy than was observed in the n case. This disagreement is a 
source of systematic error, which we try to minimize by the use of the very loose energy cut on the n shower. 




E(n) (GeV) Hits in 50-degree cone Second moment (cm 2 ) 



FIG. 8. Comparisons of distributions of selection variables for n from J/ip —¥ nn (solid line) with those from J/ip — > pnn~ 
(points): (a) deposited energy in the EMC, (b) the number of EMC hits in the 50° cone around the n shower, and (c) the 
second moment of the EMC energy deposit. 




E(n) (GeV) 

FIG. 9. Comparison of the distribution of the deposited energy in the EMC for n showers from J /if) — > nn (points) with that 
from J/ip — >• pmr + (solid line). 



For the E extra cut, we use J/ip — »• pp to obtain the efficiency (e^l*" )• The requirements are identical to those 
described in Sect. IV Al Our selection of J/ip — > pp does not depend on information from the calorimeter, so the 
behavior of pp in the EMC can be used to verify the efficiency of the E extra cut for J/ip — > nn. Figure [TU] shows the 
comparison of the E extra distributions for J/ip — > pp and J/ip —> nn. We require the angle between the n and n to 
be greater than 177° to suppress background for this comparison. We find that the proportion of E extra = events 



12 



in J / ip — > pp and J ftp nn agree well. The ratio of J /if) — > pp events with or without the requirement E extra — is 
calculated as the efficiency of the E extra cut. 




UO 0.2 0.4 0.6 0.8 1.0 



FIG. 10. Comparison of the distribution of E ex tra in J/i/) — > nn (points) with that from J /if) — » pp (solid line). 



Finally, we determine the overall efficiency from the product of the three component efficiencies described above: 



data data data 

t — t fi e n e E e:Dtra , ■ 

To facilitiate measurement of the angular distribution as well as the yield, the product efficiency is determined in 16 
bins in cos# (cosine of the n polar angle) from -0.8 to 0.8. Figure [TT] shows the efficiency for each cut and the product 
as a function of cos6>. The loss of n efficiency near cos 6* = ±0.8 is caused by the requirement on the number of EMC 
hits in a 50° cone around the shower. To smooth the bin-to-bin statistical fluctuations in the efficiency correction, we 
fit with a fifth-order polynomial (Fig. ITTT d)). 
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FIG. 11. Component selection efficiencies for J/ip — > nn as a function of cos6>: (a) n selection, (b) n selection, and (c) E ex t ra 
cut; (d) overall product efficiency as computed (points) and smoothed by fitting to a fifth-order polynomial (line). 



C. Angular Distribution and branching fraction 

The number of J/ip — > nn events in each cos 9 bin is obtained by fitting the distribution of the angle between the n 
and n. The signal shape is determined with the J/ip — > pp sample. Because of the 1.0- T magnetic field in the BESIII 
detector, the angular distribution in J/ip — s- pp must be corrected before being applied to J/ip —> nn. The signal 
shape (^) for fitting the distribution of the angle between n and n can be expressed in terms of d9 p (d8p) and d(p p 
(d<pp), the polar and azimuthal angles between the shower position and the extrapolated EMC position of p (p), as 
follows: 



* = cos 2 9 p )d(f>l p + dOlp, 

where 9 p is the polar angle of the proton track, d9 p p = d9 p + d6p and d<p p p — depp + d<pp. The background shape is 
fixed to the shape of the inclusive background, while signal and background normalizations are allowed to float in 
each cos 9 bin. A sample fit for one cos 9 bin is shown in Fig. [T2l 

After obtaining the bin-by-bin signal yields, we fit the resulting cosf? distribution with the function ^4(1 + 
a cos 2 #)e(cos 9), where A gives the overall normalization and e(cos6>) is the corrected efficiency. The resulting 
angular distribution and fit are shown in Fig. 1131 The \ 2 f° r the fit is 13 for 14 degrees of freedom, and the value 
determined for the angular-distribution parameter is a = 0.50 ± 0.04 (statistical error only). 

The raw number of J/ip nn events in the range cos 9 = [-0.8, 0.8] is N(-0.8, 0.8) = 35891 ±211. The efficiency- 
corrected yield obtained from the cos 9 fit is N cor (— 0.8, 0.8) = 354195±2078. The fitted value of a is used to determine 
the total number of J/ip — > nn, events in the full angular range of cos 6* = [—1.0, 1.0] as follows: 



f „(1 + a cos 2 , 

iV cor (-1.0, 1.0) = 7V cor (-0.8,0.8) x — jj^g— — = 466590 ± 2737. 

J_ . 8 (l + acos 2 " s 
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Combining this total yield with the number of J/ip events in our sample, we find the branching fraction to be 

B{J/ip ->• nn) = (2.07 ± 0.01) x 10~ 3 , 

where the error is only statistical. 
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FIG. 12. Fit to the angle between the n and n for cos# in [-0.3,-0.2]. 




FIG. 13. The points represent the measured distribution of cos# for the n in J/i|i -> nn candidate events, with error bars 
that are the quadratic sums of the statistical and efficiency uncertainties. The line represents the fit of the distribution to the 
functional form given in the text, and is used to determine the normalization and the angular distribution parameter a. 



D. Systematic Errors and Results 



The different n (n) momentum distributions in J/ip — > l pfVK~ (c.c.) and J/ip — > nn may introduce systematic 
uncertainties in the n (n) efficiency determination. We change the missing-momentum range from (1.1 — 1.2) GeV/c 
to (1.0 — 1.1) GcV/c when selecting the n (n) sample from J/ip — > phn~ (c.c.) and take the resulting differences 
in a and the branching fraction as systematic errors. This estimation is cross-checked with a lower- statistics sample 
obtained from a separate BESIII data sample collected at the ^(3686) resonance. High-momentum n (n) candidates 
are selected from the decay ^(3686) — > tt + tt~J/i/j, J/ip —> pnir~ (c.c). The average n (n) efficiencies obtained with 
this sample are consistent with those from J/ip — > pnn" (c.c.) within statistical errors. 

A second source of systematic error in the n (n) efficiency is the effect of the requirement that the shower be 
within 10° of the expected direction. We estimate the systematic errors due to this requirement by removing it and 
determining the changes in the results. We sum these two systematic error in quadrature to obtain the total systematic 
errors due to the selection of n (n). For the a determination, the n and n errors are 0.04 and 0.09, respectively, and 
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for the branching fraction they are 5 x 10 -5 and 1.2 x 10~ 4 . For the branching fraction the n efficiency is the largest 
source of uncertainty in our measurement. 

As we did for J/tp — > pp, we use a toy MC method to estimate systematic errors due to the statistical uncertainties in 
the efficiency. For a, this systematic error is 0.17, the largest contributor to the overall uncertainty of the measurement. 
For the branching fraction this error it is 7.1 x 1CP 5 and is the second-largest contributor. 

We change the background shape for each exclusive MC background channel and repeat the fit of the angle between 
n and n. The largest variation observed for any case considered is assigned as the systematic error. 

Our signal shape in fitting the n — h angle was obtained from J/tp — > pp. The correction of the p — p angular 
distribution into one appropriate for n — n in J/tp — > nh is a source of systematic uncertainty. To assess this we used 
a sideband subtraction instead of the fit to the angular distribution. We normalize the yield of MC background in the 
signal region (170° — 180°) by the numbers of events in the sideband range (160° — 170°) for data and MC background. 
Then we take the background-subtracted number of events in the signal region (170° — 180°) as the yield in each cos 6* 
bin. The differences between this alternative method and the standard method are assigned as systematic errors. 

The angular resolution can introduce systematic uncertainty both through the binning and through the | cos6*| < 
0.8 cut. We use the J/tp — > pp sample to evaluate the cos 9 resolution for J/tp — > nn. In the data dcos9 = 
cos 9p ext — cos 9p emc is calculated as the equivalent of the resolution in cos 9 for n, where cos 9p ext represents the 
extrapolated position of the p at the EMC and cos#p emc is the reconstructed position in the EMC. Here we assume 
that the position reconstruction of n in the EMC is similar to that for p, because both are dominated by hadronic 
interactions. To estimate the systematic error, we smear the cos 9 of n with the distribution of dcos9 and redo the 
fit in each bin, and the fit to the angular distribution. The resulting changes are taken as the systematic errors. 

For the all-neutral nn final state the trigger efficiency is another potential source of uncertainty. We correct the 
efficiency curve with the MC-determined trigger efficiency and redo the fit. The resulting changes in the branching 
fraction and a are taken as systematic uncertainties due to trigger efficiency. 

To consider the systematic error from interference between the J/tp peak and the continuum, we write the total 
cross section a nfl as 



I /„cont. i \/127rr ee r fot o i0Q\|2 
CTnn-IV^nft H 2 i ■ "H (Ai - e *J I 

s — m z + imi tot 



where E n and S are the EM and strong amplitudes for J/tp — >• nn and cp is the phase angle between them. cr^ n '' is 
the nn cross section contributed by continuum under the J/tp peak. These values are taken from Section I VIII and, 
as discussed there, the EM amplitude E n should be opposite to E p . The cross section a(e + e~ — > nn) close to J/tp is 
assumed to lie between cr(e + e~ — s- pp) and er(e + e~ — > pp) ■ (u n /u p ) 2 , so the E n , which is in proportion to y/ (T^™ 4 -, 
ranges from E p to E p (u n / u p ), where u n and u p are the magnetic moments of the neutron and proton. To estimate 
the uncertainty from the continuum, we take the larger one, c^™*' ~ cr p p nt ' , therefore also E n ~ E p . The difference 
with and without cr^™*' is assigned as systematic error. 

Finally, we change a by ±1 a (including the systematic error) and reevaluate the branching fraction to estimate 
the systematic error in the branching fraction measurement. 

Table Hill summarizes the systematic uncertainties and their sum in quadrature. The final results for our J/tp — > nn 
measurements are as follows: 



a = 0.50 ±0.04 ±0.21, and 



B(J/ip -> nn) = (2.07 ± 0.01 ± 0.17) x 10^. 

The branching fraction measurement is consistent with the previous world average @ and improves the overall 
precision by about a factor of 2.3. 



VII. SUMMARY 



We have used the world's largest sample of J/tp decays to make new measurements of the branching fractions and 
production-angle distributions for J/tp pp and J/tp — ¥ nn, obtaining the branching fractions B(J/tp — > pp) = 
(2.112 ± 0.004 ± 0.031) x 10" 3 and B(J/tp -> nn) = (2.07 ± 0.01 ± 0.17) x 10" 3 . These results represent significant 
improvements over previous measurements. The angular distributions for both decays are well described by the 
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TABLE III. Systematic errors for J/ip — > nn. 



Sources 


Effect on a 


Effect on B (10~ a ) 


n selection 


0.04 


0.05 


n selection 


u.uy 


n i o 
U.12 


Efficiency correction statistics 


n 17 


U.U f 


Background 


0.03 


0.03 


Signal shape 


0.02 


0.06 


cos 6 resolution 


0.05 


0.01 


Trigger 


0.03 


0.01 


a value 




0.01 


number of J/ip 




0.03 


continuum 




0.01 


Total 


0.21 


0.17 



functional form 1 + q cos 2 9, with measured values of a = 0.595 ±0.012 ±0.015 for J/ip — > pp, and a = 0.50 ±0.04± 0.21 
for J/ip — > nn. 

The pp angular distribution can be decomposed as \C¥^ | 2 (1 + cos 2 9) + (2M p /Mj/^) 2 \C p \ 2 sin 2 9, where C^ 1 and 
Cp are the total helicity ±1 and decay amplitudes. In terms of the angular parameter a, the ratio of amplitudes is 
\ C pl C p\ = Mj/^/{2M p )yJ{l - a) /{I + a). With our measured values for a, we find \Cf /C^\ = 0.832±0.015±0.019 
and |C„ /C~f \ = 0.95 ± 0.05 ± 0.27, respectively. These measurements permit discrimination among the different 
proposed models (2414301 ] . 

The relative phase between the strong and EM amplitudes can be obtained by comparing B(J/ip —> pp) and 
B( J/ip — > nn). The J/ip EM decay amplitudes are related to the corresponding continuum cross sections close to the 
J/ip as follows: E%(J/ip -> 7* -> NN) = B{J/ip -> fifi) ■ a(e + er -> NN)/a(e+e~ -> /i/i). Present data jsj suggest 
that cr(e + e _ — > pp) ~ (9 ± 3) pb, if fitted with a smooth W~ w as expected at high enough center-of-mass energies 
W . For cr(e + e~ — > nn) the only available data [HI, |3^| are close to threshold. In the following it is assumed that 
the neutron time-like dominant magnetic form factor is negative at these center-of-mass energies, like the magnetic 
moment, as predicted by dispersion relations [34|. The cross section a(e + e~ — > nn) close to J/ip is assumed to 
lie between a(e + e~ —> pp), as is seen in the present data close to threshold, and a(e + e~ —> pp) ■ {a n /a p ) 2 , as in 
the space- like region [351 ] . Taking into account these hypotheses and their overall uncertainties, and neglecting the 
contribution of continuum amplitudes, the strong amplitude S is given by 

S 2 = [(B(J/ip -> pp) - E 2 )E n + {B{J/iP -> nn) - E 2 )E p }/(E n + E p ) = (2.038 ± 0.094) • 10" 3 , 

and the phase cp between the strong and EM amplitudes is found to be 

<p = cos^HBiJ/ip -> pp) -S 2 - E 2 )/{2SE P )} = (88.7 ± 8.1)°. 

The uncertainty in the phase is mostly due to the B(J/ip — > nn) systematic error. This determination confirms the 
orthogonality of the strong and EM amplitudes within the precision of our measurement. 
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